Long-read sequencing reveals the structural complexity of genomic integration of HPV DNA in cervical cancer cell lines.

BACKGROUND: Cervical cancer (CC) causes more than 311,000 deaths annually worldwide. The integration of human papillomavirus (HPV) is a crucial genetic event that contributes to cervical carcinogenesis. Despite HPV DNA integration is known to disrupt the genomic architecture of both the host and viral genomes in CC, the complexity of this process remains largely unexplored. RESULTS: In this study, we conducted whole-genome sequencing (WGS) at 55-65X coverage utilizing the PacBio long-read sequencing platform in SiHa and HeLa cells, followed by comprehensive analyses of the sequence data to elucidate the complexity of HPV integration. Firstly, our results demonstrated that PacBio long-read sequencing effectively identifies HPV integration breakpoints with comparable accuracy to targeted-capture Next-generation sequencing (NGS) methods. Secondly, we constructed detailed models of complex integrated genome structures that included both the HPV genome and nearby regions of the human genome by utilizing PacBio long-read WGS. Thirdly, our sequencing results revealed the occurrence of a wide variety of genome-wide structural variations (SVs) in SiHa and HeLa cells. Additionally, our analysis further revealed a potential correlation between changes in gene expression levels and SVs on chromosome 13 in the genome of SiHa cells. CONCLUSIONS: Using PacBio long-read sequencing, we have successfully constructed complex models illustrating HPV integrated genome structures in SiHa and HeLa cells. This accomplishment serves as a compelling demonstration of the valuable capabilities of long-read sequencing in detecting and characterizing HPV genomic integration structures within human cells. Furthermore, these findings offer critical insights into the complex process of HPV16 and HPV18 integration and their potential contribution to the development of cervical cancer.

Viral integration promotes SV40T-induced immortalization by disturbing the expression of DNA/chromosome- and ECM-associated functional genes.

Lentivirus containing simian virus 40 large T antigen (SV40T) is routinely used to induce cell immortalization. However, the roles of viral integration itself in this progress is still controversial. Here, we transformed primary mouse embryonic fibroblasts (MEFs) with SV40T lentivirus and studied the roles of viral integration in the immortalization using RNA sequencing (RNA-seq) and whole genome sequencing (WGS). During the immortalization, differentially expressed genes (DGEs) are enriched in viral infection and several diverse activities. However, DEGs between immortalized and aging cells are significantly enriched in DNA/chromosome- and extracellular matrix (ECM)-associated activities. Gene regulatory network (GRN) analysis shows that although p53 is a key regulatory factor, many other transcription factors also play critical roles in the process, like STAT1. Of these DEGs, 32 genes have viral integration in their coding and/or regulatory regions. Our findings suggest that viral integration may promote SV40T-mediated immortalization by disturbing the expression of DNA/chromosome- and ECM-associated genes.

The analysis of HPV integration sites based on nanopore sequencing and the profiling changes along the course of photodynamic therapy.

OBJECTIVE: To detect the HPV genotype and integration sites in patients with high-risk HPV infection at different stages of photodynamic therapy using nanopore technology and to evaluate the treatment effect. METHODS: Four patients with HPV infection were selected and subjected to photodynamic therapy, and cervical exfoliated cell was sampled at before treatment, after three courses of treatment and six courses of treatment, their viral abundance and insertion sites were analyzed by nanopore technology, and pathological examinations were performed before and after treatment. In this study, we developed a novel assay that combined viral sequence enrichment and Nanopore sequencing for identification of HPV genotype and integration sites at once. The assay has obvious advantages over qPCR or NGS-based methods, as it has better sensitivity after viral sequences enrichment and can generate long-reads (kb to Mb) for better detection rate of structure variations, moreover, fast turn-around time for real-time viral sequencing and analysis. RESULTS: The pathological grade was reduced in all four patients after photodynamic therapy. Virus has been cleared in two cases after treatment, the virus amount reduced after treatment but not completely cleared in one case, and two type viruses were cleared and one type virus persisted after treatment in the last patient with multiple infection. Viral abundance and the number of integration sites were positively correlated. Gene enrichment analysis showed complete viral clearance in 1 patient and 3 patients required follow-up. CONCLUSION: Nanopore sequencing can effectively monitor the abundance of HPV viruses and integration sites to show the presence status of viruses, and combined with the results of gene enrichment analysis, the treatment effect can be dynamically assessed.

Prevalent and Disseminated Recombinant and Wild-Type Adeno-Associated Virus Integration in Macaques and Humans.

Integration of naturally occurring adeno-associated viruses (AAV; wild-type AAV [wtAAV]) and those used in gene therapy (recombinant AAV [rAAV]) into host genomic DNA has been documented for over two decades. Results from mouse and dog studies have raised concerns of insertional mutagenesis and clonal expansion following AAV exposure, particularly in the context of gene therapy. This study aimed to characterize the genomic location, abundance, and expansion of wtAAV and rAAV integrations in macaque and human genomes. Using an unbiased, next-generation sequencing-based approach, we identified the genome-wide integration loci in tissue samples (primarily liver) in 168 nonhuman primates (NHPs) and 85 humans naïve to rAAV exposure and 86 NHPs treated with rAAV in preclinical studies. Our results suggest that rAAV and wtAAV integrations exhibit similar, broad distribution patterns across species, with a higher frequency in genomic regions highly vulnerable to DNA damage or close to highly transcribed genes. rAAV exhibited a higher abundance of unique integration loci, whereas wtAAV integration loci were associated with greater clonal expansion. This expansive and detailed characterization of AAV integration in NHPs and humans provides key translational insights, with important implications for the safety of rAAV as a gene therapy vector.

Characterization of integrated Marek's disease virus genomes supports a model of integration by homology-directed recombination and telomere-loop-driven excision.

IMPORTANCE: Marek's disease virus (MDV) is a ubiquitous chicken pathogen that inflicts a large economic burden on the poultry industry, despite worldwide vaccination programs. MDV is only partially controlled by available vaccines, and the virus retains the ability to replicate and spread between vaccinated birds. Following an initial infection, MDV enters a latent state and integrates into host telomeres and this may be a prerequisite for malignant transformation, which is usually fatal. To understand the mechanism that underlies the dynamic relationship between integrated-latent and reactivated MDV, we have characterized integrated MDV (iMDV) genomes and their associated telomeres. This revealed a single orientation among iMDV genomes and the loss of some terminal sequences that is consistent with integration by homology-directed recombination and excision via a telomere-loop-mediated process.

Unraveling the palindromic and nonpalindromic motifs of retroviral integration site sequences by statistical mixture models.

A weak palindromic nucleotide motif is the hallmark of retroviral integration site alignments. Given that the majority of target sequences are not palindromic, the current model explains the symmetry by an overlap of the nonpalindromic motif present on one of the half-sites of the sequences. Here, we show that the implementation of multicomponent mixture models allows for different interpretations consistent with the existence of both palindromic and nonpalindromic submotifs in the sets of integration site sequences. We further show that the weak palindromic motifs result from freely combined site-specific submotifs restricted to only a few positions proximal to the site of integration. The submotifs are formed by either palindrome-forming nucleotide preference or nucleotide exclusion. Using the mixture models, we also identify HIV-1-favored palindromic sequences in Alu repeats serving as local hotspots for integration. The application of the novel statistical approach provides deeper insight into the selection of retroviral integration sites and may prove to be a valuable tool in the analysis of any type of DNA motifs.

Mitochondrial DNA is a target of HBV integration.

Hepatitis B virus (HBV) may integrate into the genome of infected cells and contribute to hepatocarcinogenesis. However, the role of HBV integration in hepatocellular carcinoma (HCC) development remains unclear. In this study, we apply a high-throughput HBV integration sequencing approach that allows sensitive identification of HBV integration sites and enumeration of integration clones. We identify 3339 HBV integration sites in paired tumour and non-tumour tissue samples from 7 patients with HCC. We detect 2107 clonally expanded integrations (1817 in tumour and 290 in non-tumour tissues), and a significant enrichment of clonal HBV integrations in mitochondrial DNA (mtDNA) preferentially occurring in the oxidative phosphorylation genes (OXPHOS) and D-loop region. We also find that HBV RNA sequences are imported into the mitochondria of hepatoma cells with the involvement of polynucleotide phosphorylase (PNPASE), and that HBV RNA might have a role in the process of HBV integration into mtDNA. Our results suggest a potential mechanism by which HBV integration may contribute to HCC development.

Oncogenic virus integration: Moving toward clinical applications.

The majority of oncogenic viruses are capable of integrating into the host genome, posing significant challenges to clinical control. Recent conceptual and technological advances, however, offer promising clinical applications. Here, we summarize the advances in our understanding of oncogenic viral integration, their clinical relevance, and the future perspectives.

HPV integration generates a cellular super-enhancer which functions as ecDNA to regulate genome-wide transcription.

Human papillomavirus (HPV) integration is a critical step in cervical cancer development; however, the oncogenic mechanism at the genome-wide transcriptional level is still poorly understood. In this study, we employed integrative analysis on multi-omics data of six HPV-positive and three HPV-negative cell lines. Through HPV integration detection, super-enhancer (SE) identification, SE-associated gene expression and extrachromosomal DNA (ecDNA) investigation, we aimed to explore the genome-wide transcriptional influence of HPV integration. We identified seven high-ranking cellular SEs generated by HPV integration in total (the HPV breakpoint-induced cellular SEs, BP-cSEs), leading to intra-chromosomal and inter-chromosomal regulation of chromosomal genes. The pathway analysis revealed that the dysregulated chromosomal genes were correlated to cancer-related pathways. Importantly, we demonstrated that BP-cSEs existed in the HPV-human hybrid ecDNAs, explaining the above transcriptional alterations. Our results suggest that HPV integration generates cellular SEs that function as ecDNA to regulate unconstrained transcription, expanding the tumorigenic mechanism of HPV integration and providing insights for developing new diagnostic and therapeutic strategies.

VIS Atlas: A Database of Virus Integration Sites in Human Genome from NGS Data to Explore Integration Patterns.

Integration of oncogenic DNA viruses into the human genome is a key step in most virus-induced carcinogenesis. Here, we constructed a virus integration site (VIS) Atlas database, an extensive collection of integration breakpoints for three most prevalent oncoviruses, human papillomavirus, hepatitis B virus, and Epstein-Barr virus based on the next-generation sequencing (NGS) data, literature, and experimental data. There are 63,179 breakpoints and 47,411 junctional sequences with full annotations deposited in the VIS Atlas database, comprising 47 virus genotypes and 17 disease types. The VIS Atlas database provides (1) a genome browser for NGS breakpoint quality check, visualization of VISs, and the local genomic context; (2) a novel platform to discover integration patterns; and (3) a statistics interface for a comprehensive investigation of genotype-specific integration features. Data collected in the VIS Atlas aid to provide insights into virus pathogenic mechanisms and the development of novel antitumor drugs. The VIS Atlas database is available at https://www.vis-atlas.tech/.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains.

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ T cells, highlighting the importance of host genome organization in HIV-1 infection.

Antiretroviral APOBEC3 cytidine deaminases alter HIV-1 provirus integration site profiles.

APOBEC3 (A3) proteins are host-encoded deoxycytidine deaminases that provide an innate immune barrier to retroviral infection, notably against HIV-1. Low levels of deamination are believed to contribute to the genetic evolution of HIV-1, while intense catalytic activity of these proteins can induce catastrophic hypermutation in proviral DNA leading to near-total HIV-1 restriction. So far, little is known about how A3 cytosine deaminases might impact HIV-1 proviral DNA integration sites in human chromosomal DNA. Using a deep sequencing approach, we analyze the influence of catalytic active and inactive APOBEC3F and APOBEC3G on HIV-1 integration site selections. Here we show that DNA editing is detected at the extremities of the long terminal repeat regions of the virus. Both catalytic active and non-catalytic A3 mutants decrease insertions into gene coding sequences and increase integration sites into SINE elements, oncogenes and transcription-silencing non-B DNA features. Our data implicates A3 as a host factor influencing HIV-1 integration site selection and also promotes what appears to be a more latent expression profile.

Intratumoral Heterogeneity and Clonal Evolution Induced by HPV Integration.

The human papillomavirus (HPV) genome is integrated into host DNA in most HPV-positive cancers, but the consequences for chromosomal integrity are unknown. Continuous long-read sequencing of oropharyngeal cancers and cancer cell lines identified a previously undescribed form of structural variation, "heterocateny," characterized by diverse, interrelated, and repetitive patterns of concatemerized virus and host DNA segments within a cancer. Unique breakpoints shared across structural variants facilitated stepwise reconstruction of their evolution from a common molecular ancestor. This analysis revealed that virus and virus-host concatemers are unstable and, upon insertion into and excision from chromosomes, facilitate capture, amplification, and recombination of host DNA and chromosomal rearrangements. Evidence of heterocateny was detected in extrachromosomal and intrachromosomal DNA. These findings indicate that heterocateny is driven by the dynamic, aberrant replication and recombination of an oncogenic DNA virus, thereby extending known consequences of HPV integration to include promotion of intratumoral heterogeneity and clonal evolution. SIGNIFICANCE: Long-read sequencing of HPV-positive cancers revealed "heterocateny," a previously unreported form of genomic structural variation characterized by heterogeneous, interrelated, and repetitive genomic rearrangements within a tumor. Heterocateny is driven by unstable concatemerized HPV genomes, which facilitate capture, rearrangement, and amplification of host DNA, and promotes intratumoral heterogeneity and clonal evolution. See related commentary by McBride and White, p. 814. This article is highlighted in the In This Issue feature, p. 799.

Human papillomavirus integration perspective in small cell cervical carcinoma.

Small cell cervical carcinoma (SCCC) is a rare but aggressive malignancy. Here, we report human papillomavirus features and genomic landscape in SCCC via high-throughput HPV captured sequencing, whole-genome sequencing, whole-transcriptome sequencing, and OncoScan microarrays. HPV18 infections and integrations are commonly detected. Besides MYC family genes (37.9%), we identify SOX (8.4%), NR4A (6.3%), ANKRD (7.4%), and CEA (3.2%) family genes as HPV-integrated hotspots. We construct the genomic local haplotype around HPV-integrated sites, and find tandem duplications and amplified HPV long control regions (LCR). We propose three prominent HPV integration patterns: duplicating oncogenes (MYCN, MYC, and NR4A2), forming fusions (FGFR3-TACC3 and ANKRD12-NDUFV2), and activating genes (MYC) via the cis-regulations of viral LCRs. Moreover, focal CNA amplification peaks harbor canonical cancer genes including the HPV-integrated hotspots within MYC family, SOX2, and others. Our findings may provide potential molecular criteria for the accurate diagnosis and efficacious therapies for this lethal disease.

Analysis of viral integration reveals new insights of oncogenic mechanism in HBV-infected intrahepatic cholangiocarcinoma and combined hepatocellular-cholangiocarcinoma.

BACKGROUND: Integration of HBV DNA into the human genome could progressively contribute to hepatocarcinogenesis. Both intrahepatic cholangiocarcinoma (ICC) and combined hepatocellular-cholangiocarcinoma (CHC) are known to be associated with HBV infection. However, the integration of HBV and mechanism of HBV-induced carcinogenesis in ICC and CHC remains unclear. METHODS: 41 patients with ICC and 20 patients with CHC were recruited in the study. We conducted HIVID analysis on these 61 samples to identify HBV integration sites in both the tumor tissues and adjacent non-tumor liver tissues. To further explore the effect of HBV integration on gene alteration, we selected paired tumors and adjacent non-tumor liver tissues from 3 ICC and 4 CHC patients for RNA-seq and WGS. RESULTS: We detected 493 HBV integration sites in ICC patients, of which 417 were from tumor samples and 76 were from non-tumor samples. And 246 HBV integration sites were detected in CHC patients, of which 156 were located in the genome of tumor samples and 90 were in non-tumor samples. Recurrent HBV integration events were detected in ICC including TERT, ZMAT4, MET, ANKFN1, PLXNB2, and in CHC like TERT, ALKBH5. Together with our established data of HBV-infected hepatocellular carcinoma, we found that HBV preferentially integrates into the specific regions which may affect the gene expression and regulation in cells and involved in carcinogenesis. We further performed genomic and transcriptomic sequencing of three ICC and four CHC patients, and found that HBV fragments could integrate near some important oncogene like TERT, causing large-scale genome variations on nearby genomic sequences, and at the same time changing the expression level of the oncogenes. CONCLUSION: Comparative analysis demonstrates numerous newly discovered mutational events in ICC and CHC resulting from HBV insertions in the host genome. Our study provides an in-depth biological and clinical insights into HBV-induced ICC and CHC.

HIV-1 Preintegration Complex Preferentially Integrates the Viral DNA into Nucleosomes Containing Trimethylated Histone 3-Lysine 36 Modification and Flanking Linker DNA.

HIV-1 DNA is preferentially integrated into chromosomal hot spots by the preintegration complex (PIC). To understand the mechanism, we measured the DNA integration activity of PICs-extracted from infected cells-and intasomes, biochemically assembled PIC substructures using a number of relevant target substrates. We observed that PIC-mediated integration into human chromatin is preferred compared to genomic DNA. Surprisingly, nucleosomes lacking histone modifications were not preferred integration compared to the analogous naked DNA. Nucleosomes containing the trimethylated histone 3 lysine 36 (H3K36me3), an epigenetic mark linked to active transcription, significantly stimulated integration, but the levels remained lower than the naked DNA. Notably, H3K36me3-modified nucleosomes with linker DNA optimally supported integration mediated by the PIC but not by the intasome. Interestingly, optimal intasome-mediated integration required the cellular cofactor LEDGF. Unexpectedly, LEDGF minimally affected PIC-mediated integration into naked DNA but blocked integration into nucleosomes. The block for the PIC-mediated integration was significantly relieved by H3K36me3 modification. Mapping the integration sites in the preferred substrates revealed that specific features of the nucleosome-bound DNA are preferred for integration, whereas integration into naked DNA was random. Finally, biochemical and genetic studies demonstrate that DNA condensation by the H1 protein dramatically reduces integration, providing further evidence that features inherent to the open chromatin are preferred for HIV-1 integration. Collectively, these results identify the optimal target substrate for HIV-1 integration, report a mechanistic link between H3K36me3 and integration preference, and importantly, reveal distinct mechanisms utilized by the PIC for integration compared to the intasomes. IMPORTANCE HIV-1 infection is dependent on integration of the viral DNA into the host chromosomes. The preintegration complex (PIC) containing the viral DNA, the virally encoded integrase (IN) enzyme, and other viral/host factors carries out HIV-1 integration. HIV-1 integration is not dependent on the target DNA sequence, and yet the viral DNA is selectively inserted into specific "hot spots" of human chromosomes. A growing body of literature indicates that structural features of the human chromatin are important for integration targeting. However, the mechanisms that guide the PIC and enable insertion of the PIC-associated viral DNA into specific hot spots of the human chromosomes are not fully understood. In this study, we describe a biochemical mechanism for the preference of the HIV-1 DNA integration into open chromatin. Furthermore, our study defines a direct role for the histone epigenetic mark H3K36me3 in HIV-1 integration preference and identify an optimal substrate for HIV-1 PIC-mediated viral DNA integration.

BET-Independent Murine Leukemia Virus Integration Is Retargeted In Vivo and Selects Distinct Genomic Elements for Lymphomagenesis.

Moloney murine leukemia virus (MLV) infects BALB/c mice and induces T-cell lymphoma in mice. Retroviral integration is mediated by the interaction of the MLV integrase (IN) with members of the bromodomain and extraterminal motif (BET) protein family (BRD2, BRD3, and BRD4). The introduction of the W390A mutation into MLV IN abolishes the BET interaction. Here, we compared the replication of W390A MLV to that of wild-type (WT) MLV in adult BALB/c mice to study the role of BET proteins in replication, integration, and tumorigenesis in vivo. Comparing WT and W390A MLV infections revealed similar viral loads in the blood, thymus, and spleen cells. Interestingly, W390A MLV integration was retargeted away from GC-enriched genomic regions. However, both WT MLV- and W390A MLV-infected mice developed T-cell lymphoma after similar latencies represented by an enlarged thymus and spleen and multiorgan tumor infiltration. Integration site sequencing from splenic tumor cells revealed clonal expansion in all WT MLV- and W390A MLV-infected mice. However, the integration profiles of W390A MLV and WT MLV differed significantly. Integrations were enriched in enhancers and promoters, but compared to the WT, W390A MLV integrated less frequently into enhancers and more frequently into oncogene bodies such as Notch1 and Ppp1r16b. We conclude that host factors direct MLV in vivo integration site selection. Although BET proteins target WT MLV integration preferentially toward enhancers and promoters, insertional lymphomagenesis can occur independently from BET, likely due to the intrinsically strong enhancer/promoter of the MLV long terminal repeat (LTR). IMPORTANCE In this study, we have shown that the in vivo replication of murine leukemia virus happens independently of BET proteins, which are key host determinants involved in retroviral integration site selection. This finding opens a new research line in the discovery of alternative viral or host factors that may complement the dominant host factor. In addition, our results show that BET-independent murine leukemia virus uncouples insertional mutagenesis from gene enhancers, although lymphomagenesis still occurs despite the lack of an interaction with BET proteins. Our findings also have implications for the engineering of BET-independent MLV-based vectors for gene therapy, which may not be a safe alternative.

HPV Integration Site Mapping: A Rapid Method of Viral Integration Site (VIS) Analysis and Visualization Using Automated Workflows in CLC Microbial Genomics.

Human papillomavirus (HPV) integration within the host genome may contribute to carcinogenesis through various disruptive mechanisms. With next-generation sequencing (NGS), identification of viral and host genomic breakpoints and chimeric sequences are now possible. However, a simple, streamlined bioinformatics workflow has been non-existent until recently. Here, we tested two new, automated workflows in CLC Microbial Genomics, i.e., Viral Hybrid Capture (VHC) Data Analysis and Viral Integration Site (VIS) Identification for software performance and efficiency. The workflows embedded with HPV and human reference genomes were used to analyze a publicly available NGS dataset derived from pre- and cancerous HPV+ cervical cytology of 21 Gabonese women. The VHC and VIS workflow median runtimes were 19 and 7 min per sample, respectively. The VIS dynamic graphical outputs included read mappings, virus-host genomic breakpoints, and virus-host integration circular plots. Key findings, including disrupted and nearby genes, were summarized in an auto-generated report. Overall, the VHC and VIS workflows proved to be a rapid and accurate means of localizing viral-host integration site(s) and identifying disrupted and neighboring human genes. Applying HPV VIS-mapping to pre- or invasive tumors will advance our understanding of viral oncogenesis and facilitate the discovery of prognostic biomarkers and therapeutic targets.

RAISING is a high-performance method for identifying random transgene integration sites.

Both natural viral infections and therapeutic interventions using viral vectors pose significant risks of malignant transformation. Monitoring for clonal expansion of infected cells is important for detecting cancer. Here we developed a novel method of tracking clonality via the detection of transgene integration sites. RAISING (Rapid Amplification of Integration Sites without Interference by Genomic DNA contamination) is a sensitive, inexpensive alternative to established methods. Its compatibility with Sanger sequencing combined with our CLOVA (Clonality Value) software is critical for those without access to expensive high throughput sequencing. We analyzed samples from 688 individuals infected with the retrovirus HTLV-1, which causes adult T-cell leukemia/lymphoma (ATL) to model our method. We defined a clonality value identifying ATL patients with 100% sensitivity and 94.8% specificity, and our longitudinal analysis also demonstrates the usefulness of ATL risk assessment. Future studies will confirm the broad applicability of our technology, especially in the emerging gene therapy sector.